Calcul du gradient de température et des variations thermiques sur dalle béton : Guide technique complet pour parking extérieur

Les infrastructures de stationnement extérieures, en particulier les dalles de parking en béton armé ou précontraint, sont soumises à des sollicitations climatiques sévères. Contrairement aux structures intérieures ou abritées, une dalle de parking extérieure subit des cycles journaliers et saisonniers de rayonnement solaire, de convection éolienne et de rayonnement thermique vers la voûte céleste.

Ces phénomènes génèrent deux types de variations thermiques majeures : une variation uniforme de la température moyenne de la dalle (responsable des dilatations et contractions globales) et un gradient de température vertical non linéaire (responsable d'effets de courbure et de contraintes de bridage internes). Pour un ingénieur en structure ou un bureau d'études en génie civil, la quantification précise de ces variations est indispensable pour prévenir la fissuration, le gauchissement (wrapping) et la ruine prématurée des joints de dilatation.

Ce guide technique détaille les hypothèses fondamentales, la méthodologie de résolution thermodynamique par l'équation de Fourier, les approches réglementaires (Eurocode 1) et propose une application numérique concrète.

---

1. Hypothèses fondamentales et modélisation physique

Pour évaluer les sollicitations thermiques, la dalle est modélisée comme un milieu continu, homogène et isotrope d'épaisseur h (m), de grande extension horizontale par rapport à son épaisseur. On adopte une approche unidimensionnelle verticale (axe z, orienté de la surface supérieure z = 0 vers la surface inférieure z = h).

Propriétés thermophysiques du béton

• Conductivité thermique (k) : Exprime la capacité du matériau à conduire la chaleur. Pour un béton de densité normale, k ≈ 1,65 à 2,5 W/(m·K).
• Masse volumique (ρ) : Généralement fixée à ρ ≈ 2400 kg/m³ pour le béton armé.
• Capacité thermique massique (Cp) : Quantité de chaleur requise pour élever de 1 K la température d'une masse de 1 kg de béton. Cp ≈ 900 à 1000 J/(kg·K).
• Diffusivité thermique (α) : Elle caractérise la vitesse de propagation de la chaleur à travers le matériau. Elle est définie par la relation :
  α = k / (ρ · Cp) [en m²/s]

Conditions aux limites (Flux thermiques de surface)

Le comportement thermique de la dalle dépend des transferts d'énergie à ses interfaces supérieure (z = 0) et inférieure (z = h).

A la surface supérieure (z = 0), le flux net q_sup est la somme de trois composantes :

1. Le flux de rayonnement solaire direct et diffus (q_s) :
   q_s = α_s · I_g(t)
   Où α_s est l'absorptivité solaire de la surface (comprise entre 0,6 pour un béton clair et 0,9 pour un enrobé ou béton balayé sombre) et I_g(t) est l'irradiance globale reçue (W/m²).
2. Le flux de convection avec l'air ambiant (q_c) :
   q_c = h_c · (T_air(t) - T(0,t))
   Où h_c est le coefficient de transfert convectif (W/m²·K), fortement dépendant de la vitesse du vent v (souvent modélisé par la formule de MacAdams : h_c = 5,7 + 3,8 · v).
3. Le flux de rayonnement infrarouge thermique (q_r) :
   q_r = ε · σ · (T_ciel⁴(t) - T(0,t)⁴)
   Où ε est l'émissivité du béton (≈ 0,9), σ est la constante de Stefan-Boltzmann (5,67 x 10⁻⁸ W/m²·K⁴), et T_ciel est la température apparente du ciel.

A la surface inférieure (z = h), le flux q_inf dépend de la nature du parking. S'il s'agit d'un parking sur terre-plein, le transfert se fait par conduction dans le sol. S'il s'agit d'un plancher de parking aérien ouvert (dalle sur poteaux), le flux est principalement convectif et radiatif avec l'environnement sous-jacent.

---

2. Méthodologie et méthodes de calcul

Le calcul d'un profil thermique s'effectue selon deux niveaux de complexité : l'approche dynamique transitoire (méthodes numériques) et l'approche réglementaire forfaitaire.

Méthode 1 : Résolution de l'équation différentielle de Fourier (Approche Transitoire)

La distribution de la température T(z,t) dans l'épaisseur de la dalle est régie par l'équation de la conduction thermique unidimensionnelle en régime transitoire, sans source de chaleur interne :

∂T(z,t) / ∂t = α · [ ∂²T(z,t) / ∂z² ]

En raison du caractère périodique (cyclique journalier) et non linéaire des conditions aux limites climatiques, cette équation ne possède pas de solution analytique simple. Sa résolution nécessite l'emploi de méthodes numériques :

• Méthode des Différences Finies (MDF) : Discrétisation de l'espace et du temps (schéma explicite ou implicite de Crank-Nicolson).
• Méthode des Éléments Finis (MEF) : Idéale pour coupler les résultats thermiques directement à un modèle de calcul de structure mécanique.

Méthode 2 : L'approche normative Eurocode 1 (NF EN 1991-1-5)

Pour s'affranchir des calculs thermodynamiques lourds, l'Eurocode 1 Partie 1-5 propose des profils de température simplifiés pour le dimensionnement structural. Le profil thermique est décomposé en :

• Une composante de température uniforme (ΔT_u) : Elle représente la variation de la température moyenne de la dalle par rapport à une température initiale de fermeture de la structure (T₀). Elle régit les variations de longueur de la dalle.
• Une composante de gradient thermique linéaire (ΔT_M) : Elle représente la différence de température entre la fibre supérieure et la fibre inférieure. Elle engendre un moment fléchissant thermique si les rotations aux appuis sont bridées.
• Une composante non linéaire (ΔT_E) : Le profil réel étant parabolique ou exponentiel, cette composante induit des contraintes d'auto-équilibrage internes (contraintes de peau).

Pour une dalle de type 1 (dalle en béton), l'Eurocode spécifie des valeurs de gradients thermiques linéaires équivalents (ΔT_heat pour l'échauffement de la face supérieure, et ΔT_cool pour le refroidissement). Pour une dalle extérieure sans revêtement isolant, les valeurs de calcul typiques de ΔT_heat oscillent entre 10°C et 15°C selon la localisation géographique.

---

3. Schéma structurel et distribution du gradient

Le schéma ci-dessous modélise la décomposition thermique d'un profil non linéaire réel en une composante uniforme (dilatation pure) et une composante de gradient linéaire (courbure thermique), conformément aux principes mécaniques de Bernoulli.

Pour illustrer le comportement d'une dalle de parking extérieure, réalisons le calcul mécanique des contraintes générées par l'apparition d'un gradient de température linéaire sur une dalle hyperstatique.

Hypothèses de l'application numérique

• Géométrie : Dalle continue coulée en place, épaisseur h = 0,20 m.
• Propriétés du béton :
  • Module d'élasticité longitudinal : E_cm = 33 GPa = 33 x 10⁹ Pa
  • Coefficient de dilatation thermique linéaire : α_t = 10 x 10⁻⁶ K⁻¹
• Sollicitation thermique (Gradient linéaire d'été) :
  • Température de la surface supérieure (z = 0) : T_sup = 42°C
  • Température de la surface inférieure (z = 0,2) : T_inf = 24°C
  • La différence de température totale est de ΔT = T_sup - T_inf = 18°C.
  • Le gradient thermique linéaire est défini par : dT/dz = ΔT / h = 18 / 0,20 = 90 °C/m.

Étape 1 : Calcul de la courbure libre théorique (χ_th)

Si la dalle était parfaitement libre de se déformer (sans liaisons ni poids propre bloquant), le gradient de température induirait une courbure uniforme, notée χ_th :

χ_th = (α_t · ΔT) / h